Gate controlled Spin-Density Wave and Chiral FFLO Superconducting phases in 

interacting Quantum Spin Hall edge states 
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We explore the phases exhibited by an interacting quantum spin Hall edge state in the presence 
of finite chemical potential (applied gate voltage) and spin imbalance (applied magnetic field). We 
find that the helical nature of the edge state gives rise to orders that are expected to be absent 
in non-chiral one-dimensional electronic systems. For repulsive interactions, the ordered state has 
an oscillatory spin texture whose ordering wavevector is controlled by the chemical potential. We 
analyze the manner in which a magnetic impurity provides signatures of such oscillations. We find 
that finite spin imbalance favors a finite current carrying groundstate that is not condensed in the 
absence of interactions and is superconducting for attractive interactions. This state is characterized 
by FFLO-type oscillations where the Cooper pairs obtain a finite center of mass momentum. 
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The quantum spin Hall (QSH)/2D time-reversal invari- 
ant topological insulator state [l|-[5| has been attracting 
much interest since its recent discovery. The bulk of this 
insulating state is gapped, and characterized by a non- 
trivial Z2 topological invariant which distinguishes this 
state from an ordinary band insulator. The most inter- 
esting low energy physics occurs at the QSH edge where 
an odd- number of counter-propagating Kramers' pairs of 
edge states exist. These edge states form a helical liq- 
uid wherein the spin polarization is correlated with the 
direction of motion. Notably, these counter-propagating 
channels are robust even when disorder is present as long 
as time-reversal symmetry is preserved. This new type of 
ID liquid can exhibit unusual features, such as, fractional 
charge, Kramers' pairs of Majorana bound states, and in- 
dividual Major ana bound states when in the proximity of 
magnets, superconductivity, or both, respectively [gNs]. 

In this work, instead of considering proximity-induced 
effects in the QSH edge, we explore the rich ID phases 
intrinsic to the helical liquid in the presence of interac- 
tions. We observe that repulsive interactions can lead 
to a spin-density wave phase that is unique to the heli- 
cal liquid and argue that it is not generated in a typical 
ID electron gas (IDEG). On the other hand, attractive 
interactions render the liquid unstable to the formation 
of a Fulde-Ferrell-Larkin-Ovchinnikov (FFLO)-type su- 
perconducting phasejol, [lo|, which is an s-wave like order 
parameter that condenses at a finite wave- vector. The 
ID nature of the helical liquid gives it an upper hand for 
hosting the FFLO-type phase, which has eluded experi- 
mental observation in higher dimensions flll-fl3l|. The two 
phases emerging from repulsive versus attractive interac- 
tions are dual to each other in that the roles of spin and 
charge are exchanged. Finally, we show that with repul- 
sive interactions a magnetic impurity acts as an effective 
experimental probe of the QSH edge, inducing oscilla- 
tions in the magnetization direction which are fundamen- 



tally different from the Ruderman-Kittel-Kasuya-Yosida 
(RKKY) oscillations of the magnetization amplitude in a 
conventional lDEG[li-[i3|. 

We begin with a heuristic analysis of the non- 
interacting QSH edge in the presence of finite chemical 
potential and spin imbalance, focusing on the fundamen- 
tal differences that give rise to new order when compared 
with a typical IDEG. As shown in Fig. [1^, the QSH edge 
consists of linearly dispersing spin-dependent modes as- 
sociated with a Dirac point centered at zero momentum, 
and is described by the Hamiltonian 



Hq = V dx 
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where v is the velocity and x is the coordinate tangent to 
the edge of the QSH sample. The operator t/jR^i^ li) (^) an- 
nihilates electron moving to the right (left) with up (down) 
spin at position x. The effects of a non-zero chemical po- 
tential and a spin imbalance can be described by 



where is an effective chemical potential for up 

(down) spin in the helical liquid. The chemical poten- 
tial, /i = ^(/i^ + /i|), can be controlled by tuning a gate 
voltage, and the spin imbalance, 5s = (i^ — ji^^ may be 
controlled by applying magnetic field in the direction per- 
pendicular to the QSH plane (or more generally, parallel 
to the direction of spin-polarization of the edge state). In 
fact, because of the spin-momentum locking on the edge, 
a spin imbalance acts to give rise to a charge current. 

Given the fundamental fields comprising the QSH edge, 
the two lowest order operators that could develop non- 
vanishing expectation values in an ordered phase are 



^l)^^^{x)^l)Li{x), Oa = tljR^{x)iljLi{x), (3) 
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These order parameters represent magnetic order {{Om)) 
and superconducting order ((Oa)) and are dual to one 
another with regards to charge and spin in that the for- 
mer carries charge and spin 2h/2 while (Om) carries 
charge 2e and spin 0. 

We now argue that for no n- vanishing /i and Ss^ these 
order parameters have the striking property that they are 
inhomogeneous in space, exhibiting oscillatory behavior 
over a characteristic length scale. We begin by tuning 
fi = 6s = and considering magnetic order. The system 
is tuned to the Dirac point and and any ferromagnetic 
order perpendicular to the spin-polarization of the edge 
states would open a gap at = since it would couple via 
a constant multiplying a Pauli spin matrix. If one tunes /i 
away from zero then, in order to open a gap at the Fermi- 
level, the magnetic order must have a finite wave- vector 
of qin^ = —2jj./v = —2kF where the superscript refers 
to the free limit (see Fig. [1^). Thus, we induce a spin- 
density wave so that a gap can open at the Fermi-level as 
opposed to a gap opening at the (buried) Dirac point for 
ferromagnetic ordering. This type of chemical potential 
driven spin-density wave is unique to the helical liquid as 
seen by noting the form of a magnetic order parameter 
for a fuh IDEG: 



X (^e^^^^V^^^(x) + e-^^^^V^^^(x)^ 



For the full IDEG the non-oscilliatory terms generically 
dominate, but these terms are completely absent for 



^Rt^Li 



non- 



the helical liquid which only has e 
vanishing. Thus, the existence of a spin-density wave 
is a unique signature of the reduced degrees of freedom 
of the helical liquid as compared to a conventional IDEG. 

Now let us consider the effects of a non-zero in the 
non-interacting limit for which we will return to the free- 
fermion Hamiltonian. In the Bogoliubov-de Gennes for- 
malism the Hamiltonian can be re- written 



BdG 



ivda:! 0(T^ - {Ss/2)r^ cr^) ^(x) 



where represents particle-hole space and spin space, 
and ^(x) = (V^i?t(^) ^Li{x) V^l^^(x) V^[^(x))^. This 
has energy levels E± = ±.\vk\ ± (55-/2 (with uncorrelated 
signs). A homogenous s-wave pairing cannot open a gap 
at the Fermi-level if ^5 7^ and is thus energetically frus- 
trated. As shown in Fig. [Dd the pairing term must have 
a finite wave- vector = 8s /v in the non-interacting 
limit to open a gap. A full IDEG would have both 

'^l?t^^^^^^^^i^^t P^^^^^S terms while the QSH edge only 
has the former. Thus, in the helical case there is always 
a ground state current of Cooper pairs in one direction 
picked by the sign of 8s since the order parameter oscil- 
lates like e*^^^ instead of cos(gA^), as in fact originally 




Figure 1. Non-interacting picture for spin density wave and 
chiral FFLO-type superconductor state formation (a) Edge 
state energy spectra at chemical potential /x. It is energeti- 
cally favorable to open a gap at the Fermi-level as opposed 
to the Dirac point, thus forming a spin-density wave with 
wavevector 2/cf. (b) Bogoliubov-de Gennes spectrum for non- 
zero 8s. Solid lines are electron states, dashed lines are hole 
states. Hybridization must occur between a solid and dashed 
line with opposite spin leading to a finite pairing wavevector 
of 8s /v. 

considered by Fulde and Ferrelljoj .We refer to this state 
as a chiral FFLO state. 

We now turn to the effects of interactions and their cru- 
cial role in determining the fate of the QSH edge state 
and energetically favorable ordered states. We derive the 
corresponding phase diagram by analyzing the form of 
the susceptibilities associated with each order and show, 
as might be expected, that magnetic (superconducting) 
order is stabilized by repulsive (attractive) interactions. 
We note that we are only considering equilibrium states. 
As in previous treatments fl7l - [l9t . we ignore Umklapp 
scattering and employ the following form for interactions 
between QSH electrons: 
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where ^2(4) represents the forward scattering strength of 
different (identical) species. These terms come directly 
from short range density-density interactions and have 
been extensively studied in Refs. [13, [l8|. As done 
previously [ITKlSl, the interacting system can be explored 
within a Luttinger liquid framework by representing the 
fermion fields in terms of boson fields (j) and 0\ iI)r-^{x) ~ 
^-i{ct^{x)-e{x)) ^ ^^^^x) - e^(^(^)+^(^)). Thus, the interact- 
ing helical liquid described hy H = Hq + + Hi is 
mapped into a free boson gas with a Hamiltonian 



dx 



K 



(5) 

where u = '^((1 + ^^)^ — (^fc)^)"^^^ renormal- 
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ized velocity and K = -, , ^gY , '^M - ) is the Luttinger 
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SDW FFLO 

(O^) = mQexp[-2inxK/u] (O^) = Aoexp[ ig^x/uK] 
\ ) 



K<1 



K=l 



K>1 



Figure 2. T — phase diagram of QSH edge for interactions 
characterized by K, finite chemical potential fi and finite spin 
imbalance 6s. 



parameter. Values of K < (>)1 represent repulsive (at- 
tractive) interactions. The chemical potential terms /i 
and Ss can be absorbed as inhomogeneous shifts of the 
bosonic fields 

4){x) = + iiKx/u, 0{x) = 0{x) - Ssx/2Ku, (6) 

which transforms the Hamiltonian to the standard form 
H = ^ j dx{uK{VeY ^ ^iV^Y). Thus, while the QSH 
system bears key differences in the physics, at the tech- 
nical level, several of its properties can be mapped to the 
extensively analyzed Luttinger liquid system describing 
the low-energy physics of a spinless interacting IDEG. 

From Eq. (|6]) , it immediately follows that the magnetic 
and superconducting orders are both associated with os- 
cillations that are renormalized by the interactions. By 
noting that Om ~ e'^'^^^^\OA ~ e^*^*^^^ and using the 
shifted forms Om ^ e-^^^^^/^O.^, Oa - e^^^^/^^^^O.^, 
we conclude that 

(Om) = mo exp [iqmx] , (Oa) = Aq exp [iqAx] (7) 

where Qm = —2iiK/u and = Ss/uK. To determine 
which of the orders dominates, we inspect the form of 
the associated susceptibilities, given by Xm/A{x^T) = 
— (T7-0mA(^, t)O^^(0, 0)), where r is imaginary time. 
We adapt the standard Luttinger liquid treatment to 
our situation to obtain the following temperature depen- 
dence in the Fourier domain: 



n2K-2 



Xm{k = Qm,^ = 0) 



(8) 



The finite wave- vector dependence reflects the oscilla- 
tory behavior in Eq. ((Tj) and the stability of a particu- 
lar order is indicated by the divergence of the associated 
susceptibility for T 0, as summarized in the phase 
diagram of Fig. [2l Hence, for repulsive interactions, 
K < 1, the system magnetically orders and is charac- 
terized by oscillations whose wave- vector Qm is controlled 
by the applied chemical potential. For attractive interac- 
tions, K > 1^ the system tends to form a superconducting 
state that shows chiral FFLO-type oscillations having the 
beautiful feature that the wave- vector ^a is completely 
tunable via an applied spin imbalance. 

Given that the currently available QSH systems are 
all in the repulsively interacting regime, we now focus 



(a) 



(b) 
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Figure 3. Magnetization oscillations around a single- magnetic 
impurity on a quantum spin Hall edge which decay as one 
moves away from the impurity. (a)/i = in the weakly re- 
pulsive regime leads to a ferromagnetic domain (b) /i > in 
the weakly repulsive regime leads to a domain with oscillatory 
magnetization direction. 



on probing the magnetic phase associated with K < 1. 
We show that a weak, localized magnetic impurity that 
provides an in-plane magnetic field H{x) = HS{x) acts 
as the simplest means of observing the oscillations in the 
magnetic order of Eq. ([7]) . The coupling to the edge liquid 
due to such a magnetic perturbation is given by 

Hh = -j^B'ipH^xHx + cryHy)S{x)ip, 

where /i^ is the Bohr magneton, ip = {ipR^ i^L\)^ 
and ^ = tan~^(i^^/i^a^). As shown in Fig. [3^, a spin 
up electron impinging the impurity effectively backscat- 
ters into a spin down electron and vice- versa. The mag- 
netic perturbation, upon suppressing the spin indices in 
the {i^R-\ i^Li) fields, exactly maps to the well-known 
quantum impurity problem in spinless quantum wires 
whose scaling properties can be easily analyzed within 
the Luttinger liquid frameworkfioj. In fact, the response 
to the impurity in our situation parallels the features 
of Friedel oscillations in the vicinity of a non-magnetic 
impurity in a spinless Luttinger |21]. At high energies 
and short distances, set by the bare magnetic impu- 
rity strength, the impurity can be treated perturbatively. 
Meanwhile at low energies and large distances, interac- 
tions renormalize its strength and the behavior is gov- 
erned by the strong coupling fixed point wherein the im- 
purity effectively splits the system into two pieces. The 
resultant magnetization in the helical liquid m+(x,t) = 
mx{x^t) -\-imy{x^t) = 2/i^(0^(x, t)) takes the form 



f{x,T,K,\li\ 



(9) 



where a is a short distance cut-off determined by the 
bulk energy gap, and / is a dimensionless decaying en- 
velope function whose form depends on the regime being 
probed jlH. For instance, in the perturbative regime, 
the susceptibility of the impurity-free system Xm deter- 
mines the response to the local impurity. For T ^ 0, this 
gives / ~ x^~'^^ for a <C X <C xo, where xq is a charac- 
teristic scale set by the bare impurity strength. On the 
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other hand, for x Xq^ the strong couphng analysis gives 



f ^ X . For a more general form of magnetic quantum 
impurity coupling, the helical liquid shows a rich range 



of behavior, including modified Kondo physics [23|, [24 1 
which, in this context, necessitates an investigation of 
the finite gate-potential induced SDW physics. 

Regardless of the strengths of the impurity and in- 
teractions, and the regimes being probed, the ubiqui- 
tous feature of the magnetization is the 2qm depen- 
dence in Eq. (j9]) that reflects SDW ordering. As il- 
lustrated in Fig. [31 the impurity thus creates oscilla- 
tions in the magnetization which decay with distance. 
More explicitly, if for example Hx ^ and Hy = we 
have {mx^rriy) ~ (cos(2^/ix), — sin(2^/ix)). To heuris- 
tically understand why this oscillation occurs, consider 
the free system {K = 1) where for a given a Fermi level 
/i, the Kramers' pair of states at the Fermi points are 
^ikpx^ t), e-^^^^l i) (Fig. [Bi). In this basis, the su- 
perpositions |±) = ^(e^^^^l t) ± e-^^^^l ;)) have 
magnetizations m+ ~ -^^-^ikpx = q. A local 

magnetic field HxS{x) breaks the time-reversal symme- 
try of the pure QSH system, removes the degeneracy be- 
tween the two states |+), |— ) since {-\-\HxS{x)a^\-\-) = 
— {—\HxS{x)a^\—) ^ 0, and forces an incomplete com- 
pensation in magnetization. We emphasize that, in con- 
trast to the oscillations that yield RKKY interactions, 
the helical nature of the QSH gives rise to spin oscilla- 
tions in direction while the magnitude remains fixed. 

In conclusion, we have shown that the helical nature 
of the QSH is unique in giving rise to ordered oscillatory 
phases in the presence of finite chemical potential and 
spin imbalance. The experimental feasibility of realizing 
and detecting the SDW phase is promising. For instance, 
in HgTe quantum wells we expect the interactions to be 
weakly repulsive and v ~ lO^m/s which leads to a char- 
acteristic tunable wavelength of around 200//i nm where 
/i is the chemical potential tuned from the edge state 
Dirac point in meV. As mentioned, depending on the 
temperature this oscillation will be modulated by a (per- 
haps strongly) decaying envelope function. To detect the 
oscillations, one could perhaps employ scanning tunnel- 
ing microscopy as has been successfully used to observe 
the RKKY oscillation [25, 26] of the magnetization near 
a magnetic impurity, and Friedel oscillations [l^l near a 
charged impurity [28]. While oscillations in the magneti- 
zation direction are harder to detect, any gate- voltage de- 
pendent oscillations would be indicative of our proposed 
SDW phase. The only modification to current setups 
would be the necessity of a back gate so that the oscilla- 
tions could be accessed. Finally, the possibility to induce 
attractive interactions in the QSH system, as has been 
achieved in ID cold atomic gases, would open up the fas- 
cinating prospect of realizing the chiral FFLO oscillatory 
superconducting phase. 
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